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Summary 


A program is underway to improve the efficiency of a three-dimensional, unsteady, Navier- 
Stokes code and generalize it for nozzle and turbopump geometries. Code modifications have 
included the implementation of parallel processing software, incorporation of new physical 
models and generalization of the multi-block capability to allow the simultaneous simulation 
of nozzles and turbopumps. The modified code has been applied to relevant nozzle and 
turbopump configurations. 

The following results are described in the 1997 Final Report: 

• Details of code modifications 

• Numerical results for a nozzle configuration 

• Numerical results for interacting nozzle/blade/vane geometries 

• Parallelization of the Navier- Stokes code using the Message Passing Interface software 
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m - Mass flow rate 

M - Mach number 

P - Static pressure 

Pt - Total pressure 
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T t - Total temperature 
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SUBSCRIPTS 
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oo - Free stream 
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- Relative frame quantity 
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Introduction 


Flow unsteadiness is a major factor in turbine performance and durability. This is es- 
pecially true if the turbine is a high work design, compact, transonic, supersonic, counter 
rotating, or uses a dense drive gas. The vast majority of modern rocket turbine designs 
fall into these categories. For example, the (Space Transportation Main Engine) STME fuel 
turbine, a high work, transonic design, w r as found to have an unsteady interrow' shock wdiich 
reduced efficiency by 2 points and increased dynamic loading by 24 percent. The Revolu- 
tionary Reusable Technology Turbopump (RRTT), which uses full flow oxygen for its drive 
gas, was found to shed vortices with such energy as to raise serious blade durability concerns. 
In both cases, the sources of the problems were uncovered (before turbopump testing) with 
the application of validated, unsteady computational fluid dynamics (CFD) to the designs. 
In the case of the RRTT and the Alternate Turbopump Development (ATD) turbines, the 
unsteady CFD codes have been used not just to identify problems, but to guide designs 
which mitigate problems due to unsteadiness. Using unsteady flow analyses as a part of the 
design process will lead to turbine designs with higher performance (which effects tempera- 
ture and mass flow rate) and fewer dynamics problems. One drawback to three-dimensional 
unsteady flow analyses, however, is that computation times are often on the order of weeks 
or months. This time scale does not allow' for parametric design studies. The objective of the 
current program is to generalize and improve the efficiency of a three-dimensional, unsteady, 
Navier-Stokes code. 


Code Improvements 

The following improvements have been made to the CORSAIR analysis during the last 
quarter: 

• A parallel version of the flow analysis has been developed using the Message Passage 

Interface (MPI) software. ° 

• The inlet boundary conditions have been modified to allow the modelling of discrete 
nozzles. This is accomplished by applying the boundary conditions on a point- wise 
basis. The inlet boundary conditions can now be specified to be solid wall, subsonic 
inflow or supersonic inflow. Presently, it is assumed that the nozzles are circular or 
elliptical in shape. It should be noted that the turbulence quantities are still being 
calculated based on the airfoil and endwall surfaces. Eventually, the solid wall .portions 
of the inlet should also be used to influence the turbulence quantities. 

• The technique for specifying the initial flow field in nozzle simulations has been mod- 
ified. For normal blade row calculations flow quantities are linearly interpolated be- 
tween the inlet and exit. In nozzle simulations (where it is assumed that the throat 
is choked) the flow quantities are interpolated between the inlet and the throat, and 
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again between the throat and the exit. Note that for choked flow the throat quantities 
are known, and the area ratio can be used to check the Mach number at the inlet. 
Application of this technique has been shown to reduce the CPU time by more than 
50% compared to initializing the flow field by interpolating the flow quantities between 
the inlet and exit. 

• A low-memory version of the code has been written in which many of the subroutines 
not necessary for nozzle simulations have been eliminated. This, in turn, has allowed 
the elimination of considerable branching logic. For a given problem size the low- 
memory version of the requires 40% less memory, and is 30% faster, than the standard 
analysis for nozzle simulations 

• PLOT3D files can now be specified to be formatted or unformatted. 

• A subroutine to mass average the flow quantities at the inlet and exit of each blade row 
has been added to the analysis to complement the existing area- averaging subroutine. 

• A comprehensive and updated version of the CORSAIR User’s Manual has been com- 
pleted. 


Numerical Simulations 
FASTRAC Nozzle Operating in RPl 

The relevant turbine nozzle geometry in this investigation is from the FASTRAC pro- 
gram [1]. Two views of the nozzle grid (constant i- and constant k-planes, respectively) 
are shown in Figs. 1 and 2. The computational grid contains 131 points in the streamwise 
(i) direction, 81 grid points in circumferential (j) direction and 21 points in the radial ( k ) 
direction. Thus, the computational grid contains a total of 222,831 grid points. The noz- 
zle is designed to have an inlet Mach number of Mi ftj 0.23 and an exit Mach number 
of M 2 ~ 2.13. The inlet static temperature is 1595° R and the ratio of specific heats is 
7 = 1.108. The Reynolds number was set at Re — 1 x 10 6 . 

Mach number contours at the inlet and exit of the nozzle are shown in Figs. 3 and 4, 
respectively. The average Mach number in the inlet plane is approximately M\ = 0.23, 
which is in excellent agreement with the design intent. The average Mach number in the 
exit plane is approximately M 2 = 2.18, which compares favorably with the design intent. 
Figures 5 thru 9 illustrate Mach contours along radial sections of the nozzle, beginning near 
the centerline and progressing towards the outer casing. These figures indicate that the 
highest Mach numbers are located near the centerline. Closer to the outer casing (Figs. 8 
and 9) the Mach number decreases due to the endwall boundary layer. Figure 10, which 
shows Mach contours in a circumferential (j = 20) plane, highlights the thickening of the 
endwall boundary layer downstream of the throat. 
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The predicted nozzle loss, which is defined as 


uj = 


Pt 


inlet 


Pt 


exit 


Pt, 


nlet 


(i) 


was u = 0.048 using the third-order accurate upwind scheme [2, 3]. It is interesting to note 
that repeating the calculation using a second-order accurate upwind scheme the predicted 
loss was u> = 0.061. Although the losses differed with order of accuracy, the average inlet/exit 
Mach numbers and flow patterns were similar. 


FASTRAC Nozzle Operating in Nitrogen 

A second simulation was performed using nitrogen instead of RPl as the operating fluid. 
Operating in nitrogen, the nozzle is designed to have an inlet Mach number of % 0.23 
and an exit Mach number of M 2 « 2.40. The inlet static temperature is 523° R and the 
ratio of specific heats is 7 = 1.467. Again, the Reynolds number was set at Re = 1 x 10 6 . 
The grid topology is identical to that used in the RPl simulation. 

Mach number contours at the inlet and exit of the nozzle are shown in Figs. 11 and 
12, respectively. The average Mach number in the inlet plane is approximately M\ = 0.22, 
which agrees closely with the design intent. The average Mach number in the exit plane is 
approximately M 2 = 2.42, which compares favorably with the design intent. Figures 13 thru 
17 illustrate Mach contours along radial sections of the nozzle, beginning near the centerline 
and progressing towards the outer casing. Qualitatively, the results are similar to those 
predicted in the RPl simulation. Figure 18, which shows Mach contours in a circumferential 
C? — 20) plane, highlights the endwall boundary layer development downstream of the throat. 

The predicted loss using the third-order accurate scheme is w = 0.071, which is higher 
than in the RPl simulation. Repeating the analysis with the second-order accurate scheme 
the losses were predicted to be u = 0.077 


FASTRAC Nozzle/Blade/Vane Geometry Operating in RPl 

A series of numerical simulations have been performed for the FASTRAC turbine con- 
figuration consisting of a modelled nozzle exit, a rotor row and a vane row. The actual 
configuration consists of 26 nozzles, 147 rotor airfoils and 75 vane airfoils. In the current 
effort, two blade count approximations have been tested: 

1. l-nozzle/6-rotors/3- vanes - Using this approximation it is assumed that there are 150 
rotor airfoils, and the rotors are scaled by 147/150 to maintain the pitch-to-chord ratio. 
The modelling of the nozzle is then exact because 6 rotor passages will completely 
encompass one nozzle. 
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Grid Type 

Inlet (Nozzle) 

Rotor 

Vane 

0 

- 

6x91x21x31 

3x91x21x31 

H 

1x14x51x31 

6x64x21x31 

3x82x31x31 

Total Points 

22,134 

605,430 

414,129 


Table 1: Grid dimensions for the l-nozzle/6-rotor/3-vane simulation. 


Grid Type 

Inlet (Nozzle) 

Rotor 

Vane 

0 

- 

2x91x21x31 

1x91x21x31 

H 

1x14x51x31 

2x64x21x31 

1x82x31x31 

Total Points 

22,134 

201,810 

138,043 


Table 2: Grid dimensions for the l-nozzle/2-rotor/l-vane simulation. 


2. l-nozzle/2-rotors/l-vane - This scaling also assumes there are 150 rotor airfoils, how- 
ever the nozzles now are approximated as a continuous slit around the annulus (see 
Fig. 19). 

The grid densities (number of passages xi xj xk) for the two simulations are presented in 
Tables 1 and 2, respectively. The total number of grid points in the l-nozzle/6-rotor/3-vane 
simulation is 1,041,693, while 361,987 grid points are being used in the l-nozzle/2-rotor/l- 
vane calculation. Figures 20 and 21 illustrate ( x — y ) and ( z — y ) views of the grid used 
in the l-nozzle/2-rotor/l-vane simulation. The average value of y + , the non-dimensional 
distance of the first grid line above the surface was approximately 1.5 for the airfoils surfaces 
and 5.0 for the endwall surfaces. A value of y + « 5.0 indicates that increased spanwise grid 
resolution may be needed. 

The turbine is designed to have an inlet Mach number of Mi w 2.11. The inlet static 
temperature is 1291° R and the ratio of specific heats is 7 = 1.108. The Reynolds number 
was set at Re = 1.4 x 10 6 . The results described below are for the l-nozzle/2-blade/l-vane 
simulation, which has completed 16 global cycles. A global cycle refers to two rotor airfoils 
passing one vane airfoil. 

Fifty snapshots of the turbine flow field were recorded during the course of one global 
cycle for the purpose of making an animation. Figures 22 thru 26 show instantaneous Mach 
number contours at 0%, 20%, 40%, 60% and 80% of a global cycle, respectively. Figures 27 
to 31 illustrate instantaneous static pressure contours at the same 5 time increments, while 
Figs. 32 thru 36 show instantaneous entropy contours. Note that the small breaks in the 
contour lines at the zonal boundaries are due to limitations in the graphics package and not 
to interpolation errors in the flow analysis. The graphics package included in the flow analysis 
first plots the contour lines on the H-grid, then “blanks” out the appropriate portion and 
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plots the O-grid contours. Thus, the H- and O-grid contours are generated independently. 
Based on the snapshots of the flow field shown in Figs. 22 to 36 the following observations 
have been made: 

• The absence of strong shock or expansion waves far upstream of the rotor suggest that 
the rotor is operating at the unique incidence condition. 

• There is a detached bow shock which stands approximately 3-4% of the rotor axial 
chord upstream of. the rotor leading edge. The shock extends into the suction side of the 
rotor passage w’here it interacts with airfoil the boundary layer. The shock/boundary 
layer interaction induces flow separation at approximately 15% of the axial chord on 
the rotor suction surface. 

• The separated flow region on the rotor suction surface extends axially to the trailing 
edge of the airfoil and circumferentially to almost mid-passage (also see Fig. 37). The 
flow within the separated region is subsonic, but the flow in the remainder of the 
passage is supersonic. The entropy contours indicate there is a significant amount of 
loss generation associated with the separated flow region. 

• There is minimal unsteady interaction between the rotor and vane airfoils because of 
the relatively large axial gap between them. 

• The flow along the suction surface of the vane stays attached until immediately up- 
stream of the trailing edge. There is a shallow separation bubble originating from the 
leading edge of the pressure surface and extending to approximately 30% of the axial 
chord (also see Fig. 38). 

• There is a normal shock at approximately 53% of the axial chord on the suction surface 
of the vane. The shock extends a maximum of about 25% of the distance across the 
vane passage; the flow does not choke at any time during the course of a global cycle. 

In an effort to understand the development of the losses within turbine, instantaneous 
total pressure contours were plotted in axial planes 30% of the rotor axial chord upstream 
of the rotor passage, midway between the rotor and vane passages, and 47% of the vane 
axial chord downstream of the vane (see Figs. 39 to 41). Note that the free stream value 
of the non-dimensional total pressure is Pt/P , ^ = 9.12 (see Table 3). Upstream of the 
rotor the contours indicate the presence of relatively thick endwall boundary layers (see 
Fig. 39). The reason for the thickness of the boundary layers is probably two-fold: 1) the 
solid wall portions of the nozzle at the inlet are located adjacent to the endwalls, and 2) the 
spanwise grid density needs to be increased. To address the problem of spanwise grid density 
requirements, a simulation has been initiated using the same streamwise and circumferential 
grid point distributions, but with 51 spanwise planes instead of 31 spanwise planes. Moving 
downstream to the region between the rotor and vane it is observed that the total pressure 
is only about 10-20% of its upstream value (see Fig. 40). There are two regions of low total 
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Variable 

INLET 

ROTOR 

VANE 

M tn 

2.11 

1.93 

0.65 

M in 

2.11 

1.38 

0.65 

M 0 ut 

1.93 

0.65 

0.39 

Mout 

1.93 

0.71 

0.39 

Tt 

1601° R 

1598° R 

1453° R 

Tt in 

1601° R 

1469° R 

1453° R 

Tt out 

1598° R 

1453° R 

1447° R 

Tt out 

1598° R 

1468° R 

1447° R 

(«/a o),» 

9.12 

7.91 

1.51 

(«/a ok 

9.12 

3.20 

1.51 

(«/*»)_ 

7.91 

1.51 

1.42 

(Pt/P~ L„ 

7.91 

1.74 

1.42 

71 

72.0° 

73.7° 

'-8.3° 

Pin 

72.0° 

69.0° 

-8.3° 

&out 

73.7° 

-8.3° 

-6.6° 

Pout 

73.7° 

-36.0° 

-6.6° 

4*i n 

0.0° 

-8.6° 

-3.2° 

i’in 

0.0° 

-5.0° 

-3.2° 

4*out 

-8.6° 

-3.2° 

-1.3° 

$ out 

-8.6° 

-0.8° 

-1.3° 

rh ( lb m jsec ) 

7.95 

7.98 

7.90 

W ( BTU/lb m ) 

- 

89.4 

- 

Vtt 

- 

0.575 

- 

Its 

- 

0.508 

- 


Table 3. Area-averaged flow quantities - high flow. 


pressure associated with the rotor wakes, but also two regions near the shroud endwall where 
the total pressure is significantly higher. The source of these high total pressure regions is 
currently under investigation. Downstream of the vane there are loss regions associated with 
the vane wake and the endwall boundary layers (see Fig. 41). 

. . Tabl ® 3 contains the circumferentially and radially area-averaged flow variables at the 
inlet and exit of each blade row. Note that at the inlet only the flow coming through the 
nozzle, and not the solid wall values, are included in the average. Based on the interrogation 
of the solution, it was determined that the mass flow rate had to be reduced to match the 
intended design-flow conditions. The mass-flow was too high in the first simulation because 
the limited spanwise grid density resulted in a modelled nozzle that was slightly too large. 
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Variable 

INLET 

ROTOR 

VANE 

M in 

2.11 

1.90 

0.45 

M tn 

2.11 

1.33 

0.45 

M ou t 

1.90 

0.45 

0.39 

Mout 

1.90 

0.78 

0.39 

Tt in 

1601 0 R 

1595° R 

1434° R 

ebuhhh 

1601° R 

1470° R 

1434° R 

Ti 

± L out 

1595° R 

1434° R 

1430° R 

Ttout 

1595° R 

1465° R 

1430° R 

(Pt/P ooL 

9.12 

7.40 

1.45 

EmilSlHH 

9.12 

3.00 

1.45 

(Pt/P ~) out 

7.40 

1.45 

1.42 

BB&iSm 

7.40 

1.82 

1.42 

^171 

72.0° 

73.8° 

-18.5° 

Pin 

72.0° 

69.4° 

-18.5° 

®out 

73.8° 

-18.5° 

-8.6° 

flout 

73.8° 

-61.0° 

o 

CO 

1 

Ain 

0.0° 

i 

3.7° 

1.3° 

A in 

0.0° 

-6.0° 

1.3° 

flout 

3.7° 

1.3° 

-1.2° 

flout 

3.7° 

6.2° 

-1.2° 

m (lb m / sec) 

7.09 

7.01 

7.16 

IV ( BTUjlbm ) 

- 

104.2 

- 

Vtt 

- 

0.644 

- 

T}ts 

- 

0.612 

- 


Table 4: Area-averaged flow quantities - design flow. 


Therefore, a simulation was performed in which the nozzle semi-minor axis was adjusted 
so that the mass flow accurately reflected the design intent. Table 4 contains the circumfer- 
entially and radially area-averaged flow variables at the inlet and exit of each blade row at 
the design-flow conditions. 

FASTRAC Nozzle/Blade/Vane Geometry Operating in Nitrogen 

A l-nozzle/2-rotor/l-vane simulation of the turbine (described above) using nitrogen as 
the operating fluid has been initiated. Th solution has completed one global cycle. 
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Implementation of Parallel Processing Software 

A parallel version of the flow analysis has been developed using the Message Passing 
Interface (MPI) software. The following logic has been utilized in the parallel analysis: 

• A blade passage is defined as the combination of the O-grid, H-grid and clearance grid 
(if applicable) used to discretize a region of the flow. 

• Depending on the number of processors available the problem is decomposed in the 
following manner: 

- if the number of processors available is greater than the number of blade passages, 
then the problem is decomposed on a passage-bv-passage basis 

- if the number of processors available is greater than the number of blade rows, 
but less than the number of blade passages, then the problem is decomposed on 
a row-by-row basis 

- if the number of processors is less than the number of blade rows, then the problem 
is decomposed into blocks of rows 

• The grid and restart files are read on each processor. 

• Each processor solves its given portion of the problem. 

• After each iteration the processors exchange boundary condition information , 

• After each processor completes its given portion, the appropriate information is com- 
municated back to the master process, 

• The grid and restart files are reconstructed, and the data is post-processed. 

Parallel simulations have been performed for two test cases: 1) a coarse-grid 2-vane/l- 
rotor/l-stator geometry in which the grid point densities in the different passages were 
well-balanced (Case 1), and 2) the l-nozzle/2-rotor/l-vane geometry described above (Case 
2), in which the grid point densities were not optimally balanced. The predicted results in 
for both cases were compared to the results obtained with the serial version of the analysis 
to verify the accuracy of the solutions. Table 5 contains initial timings obtained using the 
parallel version of the flow analysis. 


Conclusions 

A three-dimensional, unsteady, Navier-Stokes code has been generalized for nozzle and 
turbopump geometries. Code modifications have included the implementation of parallel 
processing software, the incorporation of new physical models and generalization of the 
multi-block capability. The modified code has been applied to relevant nozzle and turbopump 
configurations. 
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Case 

Number of Processors 

Speed-up 

1 

4 

3.27 

2 

3 

1.70 


Table 5: Initial CPU times using the parallel version of the code. 
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Figure 2. Constant k-plane view of the FASTRAC turbine nozzle. 
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Figure 4: Mach contours in the exit (i=131) plane of the FASTRAC turbine nozzle - RPl. 
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Figure 5: Mach contours in the k=4 plane of the FASTRAC turbine nozzle - RPl. 
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Figure 8. Mach contours in the k=16 plane of the FASTRAC turbine nozzle - RPl. 
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Figure 10: Mach contours in the j=20 plane of the FASTRAC turbine nozzle - RPl. 


25 



26 




Z(in) 

Figure 12: Mach contours in the exit (i=131) plane of the FASTRAC turbine nozzle - A 2 . 
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Figure 13: Mach contours in the k=4 plane of the FASTRAC turbine nozzle - N 2 . 
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Figure 18: Mach contours in the j =20 plane of the FASTRAC turbine nozzle - N%. 


33 



0.30 


o.io H 


Y(in ) 


-o.io H 


- 0.30 — 

4.40 


HUB 


SOLID 

WALL 



FLOW 


n i 1 1 r 

4.60 4.80 5.00 

Z(in ) 


5.20 


Figure 19: Inlet nozzle definition. 
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Figure 21: Grid topology for FASTRAC turbine (z-y). 
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Figure 24: Mach contours at midspan - 40% cycle. 
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Figure 25. Mach contours at midspan - 60% cycle. 
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Figure 27: Static pressure contours at midspan - 0% cycle. 
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Figure 28: Static pressure contours at midspan - 20% cycle. 
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Figure 29: Static pressure contours at midspan - 40% cycle. 
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Figure 30: Static pressure contours at midspan - 60% cycle 
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Figure 31: Static pressure contours at midspan - 80% cycle. 
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Figure 32: Entropy contours at midspan - 0% cycle. 
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Figure 33: Entropy contours at midspan - 20% cycle. 
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Figure 34: Entropy contours at midspan - 40% cycle. 
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Figure 35: Entropy contours at midspan - 60% cycle. 
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Figure 36: Entropy contours at midspan - 80% cycle. 
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Figure 37: Relative-frame velocity vectors in the rotor - 80% cycle. 
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Figure 38: Relative-frame velocity vectors in the vane - 80% cycle. 
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Figure 41: Total pressure contours 47% of the vane axial chord downstream of the vane 
80% cycle. 




